Matter evolution in Burgulence 



Ilya a. Bogaevsky 



O 

o 

(N 



(N 



> 
m 

o 
l> 
o 

o 

J: ■ 

■ 



X 



Independent University of Moscow 
Bolshoi Vlas'evskii per. 11, Moscow 119002, Russia 
E-mail: bogaevsk@mccme.ru 



Abstract 

In inviscid solutions of the forced Burgers equation the 
matter accumulates in the shock discontinuities. We 
describe the limit motion of particles everywhere in- 
cluding the shocks as the trajectories of a discontinu- 
ous velocity field being a generalization of the gradient 
of the limit potential. The latter is not differ entiable 
but satisfies some convexity properties which guarantee 
the existence of the gradient. It turns out that for such 
discontinuous gradient ordinary differential equations 
there are natural existence, uniqueness, and continuity 
theorems. These general results are applied for inves- 
tigation of formation and motion in plane of massive 
points which are interpreted as various clusters in the 
adhesion model of the Universe. 

Keywords: Burgers equation, shocks, massive points, 
clusters, singularities, transitions. 



1 Introduction 

The subject of this paper is the matter evolution in 
limit potential solutions of the Burgers equation with 
vanishing viscosity and external potential force. The 
Burgers equation is just the Navier-Stokes equation 
without the pressure term - its theory is well described 
in the survey 

It is well known that in such inviscid potential so- 
lutions there can be shocks, i.e. velocity discontinu- 
ities. They appear even if the initial condition and 
the external force are smooth - here and further this 
term means infinite differentiability. Generically in this 
smooth case shocks are smooth hypersurfaces with pre- 
scribed singularities. "Generically" means that other 
singularities can be killed by arbitrarily small pertur- 
bation of the smooth initial condition. 



In plane such a generic shock is a smooth curve with 
triple nodes and end points looking like shown in Fig- 
ure^ It can experience the transitions which are shown 
in Figure|21 
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Figure 1. Example of shock in plane 

Limit potential solutions of the Burgers equation de- 
scribe motion when the particles cannot pass through 
each other and adhere on shocks. It happens because a 
particle cannot leave the shock. More precisely, a par- 
ticle trajectory ending outside of the shock hes outside 
of the shock as well, but a trajectory beginning outside 
of the shock can end on the shock. In other words, the 
matter accumulates in the shocks where the density is 
infinite. This is the so-called adhesion model of mat- 
ter evolution in the Universe describing the formation 
of cellular structure of the matter (see, for example, 
PI) - the adhesion of particles is a result of interaction 
between them described by the vanishing viscosity. 

But what is the motion of particles on shocks? This 
question is answered by the present paper. Its main 
results are the following. 

1) When the viscosity is positive the trajectory of 
any particle is well defined because the velocity field is 
smooth. It turns out there exists a limit of the trajec- 
tory as the viscosity vanishes. Such limit trajectories 
describe the motion of particles in the inviscid solu- 
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Figure 2. Generic transitions of shocks in plane 



tion. The following uniqueness theorem is true: there 
is only one limit trajectory beginning at a given point 
but there can be a few limit trajectories ending at it. 
Of course, a few trajectories can end only at a point of 
the shock. 

2) How to find the limit trajectory of a given parti- 
cle? It turns out that the limit trajectory is a solution 
of the Cauchy problem for a velocity field defined by 
the limit potential. The velocity field is discontinuous 
but, nevertheless, the Cauchy problem has a unique 
solution. 

Besides, the above general results are applied for 
investigation of formation and motion in plane of mas- 
sive points (points with positive mass) which are inter- 
preted as various clusters in the adhesion model of the 
Universe. It is natural to assume that on the shock 
outside of its singularities the matter is distributed 
with positive linear density and the nodes are massive 
points. However, this idea is not correct! 

Indeed, the velocity field of the limit trajectories is 
smooth on the shock outside of its singularities and a 
cluster cannot appear there. But a cluster cannot ap- 
pear at a node with an acute angle too because the 
field of the relative (with respect to the node) veloci- 
ties looks like shown in Figure|31on the left. So particles 
pass through such a node and the matter does not ac- 
cumulate here. Otherwise, if all angles of a node are 
obtuse then the matter is trapped at the node and it is 
a growing cluster as it is shown in Figure|3|in the mid- 
dle. Its right side shows that a cluster cannot appear 
at an end point of the shock as well. 

So, a cluster is born when an acute node turns into 
an obtuse one - see the left transition in Figure^J After 
the opposite transformation the cluster stops growing 
and leaves the node - this is the left transition in Fig- 
ure|3 (In our figures such stable clusters are shown by 
white disks, growing clusters - by black disks.) A sta- 
ble cluster travels along the shock and, in particular. 




^' Jc'' 

— < — 



Figure 3. Velocities of particles around acute node, 
growing cluster at obtuse node, and end point 



can pass through an acute node and be absorbed by a 
growing cluster (the left transitions in Figures El and [7| 
respectively) . 




Figure 4. Cluster is born and starts growing 



The transitions in the middle and on the right from 
Figures 0HZI show what happens with clusters when 
transitions of shocks occur. Clusters can be involved 
only in the fifth and sixth transitions from FigureEl be- 
cause all nodes of the other ones are acute. Besides, 
generically a stable cluster cannot come to a transition 
but can appear after it as shown in Figure [S] and [T] 
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Figure 5. Cluster stops growing and leaves node 



Figure 7. Two clusters adhere 
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Figure 6. Cluster travels through nodes 



I am very grateful to U. Frisch for calling my atten- 
tion to the problem as well as to him, J. Bee, M. Blank, 
K. Khanin, R. Mohayaee, and A. Sobolevsky for fruitful 
discussions. 

2 Existence of limit trajectories 

So, we consider a material d-dimensional medium 
whose velocity is potential and described by the Burg- 
ers equation with the potential force term: 

vj' -I- (v"^ • V) v'' = -VC/ + vIS.-w'' 

r(x,0) =^o(x) 

where x G M"* is a point of the medium, v''(x, i) is 
the velocity at the point x at the time > is 
the viscosity of the medium, V = • • • , dx^ is the 
usual V-operator in M'*, and A = V- V is the Laplacian. 
The potential tA'^ of the velocity field v"^ is defined with 
respect to a function of time which can be chosen so 
that the following equation is satisfied: 



The force potential \J and the initial condition (^o 
are assumed to be smooth. (Everywhere in the present 
paper it means infinite differentiability.) Let if be the 
limit solution as the viscosity vanishes. 



(^(x, t) 



lim t). 



(2) 



(1) 



According to Theory of PDE, the potential ip'^ is 
smooth if > and t > Q. As it has been men- 
tioned before, the limit potential is continuous, but its 
gradient field can have discontinuities (shocks). 

We consider the periodical case. It means that all 
data - the given force potential C/, the initial condition 
lySOi the velocity field v"^, its potential ■0'^, and the limit 
potential - are assumed to be space-periodical. In other 
words, we consider our equations on torus. 

Remark. The periodicity requirement is techni- 
cal. Informally speaking, it is needed to guarantee that 
nothing goes to the infinity and nothing comes from the 
infinity for a finite period. Otherwise, let this informal 
requirement be satisfied and we want to apply our re- 
sults in a finite domain of the space-time. Then we can 
always consider our data as periodical with an enough 
large period. Indeed, if we change them enough far we 
will not be able to observe the influence in our flnite 
domain. 

Let : [0, +oo) be the trajectory beginning 

at an initial point a G R'^, that is the solution of the 
following Cauchy problem: 

r{t)^yr{y''{t),t), 2;''(0) = a 

which has a unique solution because the right side of 
the ordinary differential equation is a smooth space- 
periodical vector field. 

Theorem 1. Existence: For any initial point a 
there exists a limit trajectory x : [0, +oo) R'' 

x(0 = hmy^(i), x(0) = y'^(O) = a. 
1^— »o 
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The convergence is uniform on any segment [0, T] and 
the limit trajectory is continuous. 

Uniqueness: // two limit trajectories pass through 
the same point at some time then they coincide after 
that time: 

Xl{t^,) — X2{t^,) Xi{t) = X2{t) V<>i*. 

(But they may not coincide before the time: Xi{t) ^ 
X2{t) for t < t^,.) 

Continuity: The point x{t) is a continuous func- 
tion of the time t and the initial point a. 

Corollary. For any point and time > 
there is a limit trajectory x : [0, +oo) — > passing 
through the point at the time: x{t<,) — x*. 

Proof. According to Theory of ODE, for any v 
there exists a trajectory y'^ such that y'^{t^:) — x^ (be- 
cause v'' is smooth). Taking into account that our 
torus is compact we can choose a sequence Vn as 
n ^ oo such that the sequence y"" (0) converges to a 
point a,t as n — > oo. For the Hmit trajectory x with the 
initial point x{0) — a» we get a:(t*) = x,. □ 

3 Differential equation for limit trajec- 
tories 

Theorem 2. The derivative 

of the limit potential along any space direction q exists 
and can be presented as the minimum of linear func- 
tions: 

</'x.t(q:0) = min p-q (3) 

where k^^t is a compact set of momenta which depends 
on the point (x, t) of the space-time. 

Besides, any limit trajectory from TheorerrJ^ satis- 
fies the differential equation 

x+{t) = u{x{t),t) 

where u(x, t) is the center of the minimal ball contain- 
ing the set kx,t and the left side is the one-way deriva- 
tive 

^ lim + 

A^+O A 

Remark. The set k^.t consists of the hmit velocities 
at the time t at points which are outside of the shock 
and tend to the point x. 



Theorem|21is proved in Sectional and can be briefly 
explained in the following way. When ly ^ the veloc- 
ity v'^(x, t) ~ V'0'^(x, t) is the solution of the following 
minimum problem: 

|q|V2-^'':,.,(q,l)^min. 

q 

It turns out that this minimum principle remains valid 
for the limit potential. Namely, the limit velocity 
u(x, t) is the solution of the same minimum problem 
for the limit potential: 

|q|V2-^;,.,(q,l)^min. (4) 
q 

Proving this principle in Sectional we do not use that 
the potentials ip'^ are solutions of the Burgers equation 
- it is only important that their second derivatives are 
uniformly bounded above. 

Remark. The principle 10} is not variational be- 
cause there is no an integral functional to be minimized 
by this principle. In fact, it just generalizes the notion 
of a gradient for a some class of non-smooth functions 
of X and t (see Sectional for details). But such the 
generalized gradient depends on the behavior of the 
function at closed times after t. It does not happen if 
the function is smooth - then we get the usual gradient 
defined completely by the first space derivatives of the 
function. 

Let us show how the principle implies TheoremEl 
It is well known that in the case when the initial condi- 
tion ipo is smooth there is the so-called minimum rep- 
resentation 

V5(x,t) -mm{F(e,x,i)} 

where _F is a family of smooth solutions of the 
Hamilton-Jacobi equation: 

Ft{^, X, t) + I VxF(e, X, i)| V2 + C/(x, t) = 

depending smoothly on a parameter ^. (In other words, 
f is a smooth function of its variables x, and t.) 
It immediately implies that the derivative of the limit 
potential along a direction of the space-time 

(p{x + Xci,t + \t) - ip{x,t) 
V^.m^V = hm 

A — >+u A 

can be presented as the minimum of some solutions 
of the Hamilton-Jacobi equation freezed at the point 
(x,i): 

<,(q,T)= min {p • q - t|p| 2/2 - C/(x, t) r} (5) 
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where Ky^ t is a compact set of momenta which depends 
on the point (x, t) of the space-time. Substituting here 
T = and comparing with the formula Q we get that 
the sets fcx,t and -fCx,* has the same convex huh. 
Applying the principle we get that the value 

|q|V2-<,(q,l) = 

= |q|V2- min {p ■ q - |p| 2/2 - [/(x, i)} = 

= max { |p - q|V2} - C/(x,t) 

attains its minimum at the center of the minimal ball 
containing the set fcx,t because the latter has the same 
convex hull as the set -fCx,*- 

Remark. The equivalence of Theorem|21 and the 
principle Q for the limit solutions of the Burgers equa- 
tion has been observed independently on the author of 
the present paper by K. Khanin and A. Sobolevsky. 

4 Proof of Theorem[T] 

Lemma 1. For any T > there exists a constant 
C(T) bounding for all i> > and < t < T the second 
derivative of the potential in any direction of the space- 
time 

Vqq < C{T), 

where Q = (q, r), q = (qi, . . .,qd), |qp + = 1. 

Proof. This is the standard maximum principle for 
the second derivative V'qq that satisfies the equation 

^QQt + ^Vq ■ ^^Q + ■ ^rQQ + UqQ = l^A^QQ 

which is a consequence of ^ and implies the inequality 

^QQt + ■ ^Vqq + B< vM^qQ 

where 

B{T) = min[/QQ(x,i), |Q| = 1, x S R'^, < t < T. 

Then the following inequality holds: 

rit + Vijf ■ Vri < lyArj, r] = ^Pqq + Bt. 

So, when the function 77 attains its maximal value for 
t £ [0, r] we get the inequality 77* < which shows 
that the maximal value can be attained only if i = 0. 
Therefore, 

rj — Bt < max {77|f=o} ^ Bt 

that means 

V'qq < niax{V'QQ|t=o} - Bt 



but max {■i/^QQ I t^o} is defined by the initial condition 
ipo and the force potential U because 

where according to 

V't" = i^AiP'' - ^ViIj" ■ Vifj" - U. 
If, for example, 

C(T)=max{^QQ} + |i?(T)|T, 

we get the statement. □ 

Lemma 2. For any T > the convergence ^ is 
uniform on R'^ x [0, T] . 

Proof. According to Lemma^the functions 
C(r)(|x|2+<2)/2„^-(x,i) 

are convex on M"^ x [0, T]. Hence, their convergence as 
v is uniform on any compact subset according to 
Theorem^ from Sectional But the functions ^'"(x, t) 
are space-periodical and their convergence is uniform 
on R'^ X [0,r]. □ 

We are proving Theorem^ Let a,e > 0. Lemmal^l 
implies that for any sufficiently small 1^ and 

|?/'(x) - 7A*(x)| < e 

where V(x) = V'*(x) = ■0''*(x,t), x e R'^, and 

t S [0, T]. We want to get a uniform upper bound for 
the square of the distance between the corresponding 
trajectories: 

R{t) = \y{t) - y,{t)\^ 
where yit) = y''{t), y,{t) ^ y"- {t), 

?/(0) = a, y,(0)=a,, |a-a,|<a. 

According to Lemma^ 

ip*{y) - ip*{y*) < VV'*(2/*) ■{y-y*) + C\y- y,p/2, 

^iy*) - Hy) < V?A(y) • (y, ~ 2/) + C |y, - 2/1 2/2. 
Adding the inequalities we get: 

-2e < -(V7A(2/) - s^My*)) ■{y-y*) + c\y~ y,|2, 

or 

(VV'(y)-VV',(y*))-(y-y*) <2£ + C|2/-y,|2, (6) 
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that gives 

i?(t) < 4e + 2Ci?(t). 

Solving the differential inequality and talcing into ac- 
count that i?(0) < a, we get: 

p2Ct -1 

R{t)<2e ^ +ae'^\ 

o 

or in the special case C — 0: 

R[t) < Aet + a. 

The inequalities give the required uniform upper 
bound. 

5 Convex functions 

Let M be a convex subset of an m-dimensional affine 
space. (It means that for any two points of M the 
segment connecting them belongs to M as well.) A 
function / : Af — > R is called convex if it satisfies the 
inequality 

f{aX + pY)<af{X)+pf{Y) 

for all a,(3 > such that a + f3 = 1 and any points 
X,Y G M. A smooth function is convex if and only 
if its second derivative along any direction is non- 
negative. 

It is well known - see, for example - that convex 
functions have many good properties: 

Theorem 3. Let M G K™ be an open convex subset 
and f : M R be a convex function. 

Continuity: / is continuous on M . 

Differentiability: / has a finite derivative along 
any direction Q at any point X £ M : 

/^(g)^^lnn^ A • 

Moreover, 

f{X + Q) = f{X) + fxiQ) + o{\Q\), Q^O 

where the derivative fxiQ) o convex homogeneous 
function of Q: 

/^(AQ) = A/^(g), A>0 

and can be presented as 

where the set 2?x(/) is convex and consists of the sub- 
differentials of the function f at the point X . 



Sub- differential boundedness: The sub- 
differentials Vxif) is uniformly bounded if X belongs 
to a compact subset of M . 

This theorem is proved in [3] (see theorems 10.1 and 
23.1). 

Theorem 4. Let M e R™ be an open convex sub- 
set and : M ^ R be a family of convex functions 
depending on a parameter v . 

Uniform convergence: // the family converges 
/.(X)= limr(X) 

then the limit function /* is convex on M and the con- 
vergence is uniform on any compact subset of M . 

Uniform derivative boundedness: // the fam- 
ily f^ is uniformly bounded on M then the family of 
the sub-differentials 'Dxif'^) is uniformly bounded if X 
belongs to a compact subset of M . 

6 Gradient differential equations 

Let M C M™ be an open convex subset; a poten- 
tial : M ^ M be the difference of a semi-definite 
quadratic form e and a convex function /: 

^{X)^e{X)~ f{X), XeM; 

T = R™ be the tangent space to M or the space of 
velocities; T* = R™* be the cotangent space to M or 
the space of momenta; and a Hamiltonian ft, : T* ^ R 
be a smooth convex function of momenta. 

The derivative of ip can be written in the following 
form: 

^'xiQ) ^ min P-Q, PeT*, QeT 

where Vx — T>x{ip) C T* is a compact convex set of 
momenta - see Theorem|21 for details. The momenta 
from the set Vx are called sub- differentials of the po- 
tential If at the point X. 

Let the convex Hamiltonian h attain its minimal 
value on "Dx at a point Px G of the convex set 
of the sub-differentials of the potential (p. 

Definition. 1) Hamiltonian form: The velocity 
Qx = h'p^ e T is denoted by Vh'fiiX) and called the 
h-gradient of the potential tp at the point X. Here h'p^ 
is the differential of the Hamiltonian h at the point Px ■ 

2) Lagrangian form: The /i-gradient Qx is the min- 
imum point of the function 

m-^'xiQ), ;(Q) =max {P.Q-/i(P)} 
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where I is the Lagrangian being the Legendre transfor- 
mation of the Hamiltonian h. 

In order to show that the Hamiltonian and La- 
grangian forms are equivalent, let us note that 

^'xiQx)^Px-Qx (7) 
where Qx — '^hf{X). Indeed, 

min P ■ Qx = Px ■ Qx 

PeVx 

because Px is a minimum point of the smooth function 
h on the convex set Vx and for any P G Vx we get 
{P -Px)- h'p^ > 0. Besides, 

h{P) >P-Q-l{Q) 

because h{P) = maxg {P ■ Q — 1{Q)}. After minimiz- 
ing we get 

h{Px) > mm P-Q- 1{Q) = (^^(Q) - 1{Q). 

But h{Px) = Px ■ Qx — l{Qx) because Qx = h'j^ that 
gives 

hiPx) = ^'x [Qx) - l{Qx) 

together with ITJ. Comparing the last two displayed 
formulas we get 

v'x {Qx)-i{Qx)>v'x{Q)-m 

that proves the equivalence of our Hamiltonian and La- 
grangian forms of defining the /i-gradient. 

Remark. If the point Px is not defined uniquely 
the /i-gradicnt docs not depend on it. Indeed, if a 
smooth convex function has the same value on a seg- 
ment, its differentials coincide at the points of the seg- 
ment. 

Theorem 5. Let the Hamiltonian he the sum of a 
linear form and a positive semi-definite quadratic form: 

h{P) = hi{P) + h2{P); 

e be a positive semi-definite quadratic form on M™ 
and if : — > M 6e a potential with bounded sub- 
differentials such that the difference e — ip is a convex 
function on M™ . 

Then the Cauchy problem 

X+{t) = Vh{^){X{t)), XiO)^Xo, 

where the left side of the differential equation is the 
one-way derivative 



has a unique global solution which depends continuously 
on the initial point Xq and the potential Lp provided that 
the quadratic form e is fixed. 

More precisely, it means that there exists a unique 
trajectory X : [0, -t-oo) satisfying the differential 

equation for any t > and the initial point X{0) = 
Xq. Moreover, if the quadratic form e is fixed then 
the point X{t) depends continuously on the time t, the 
initial point Xq, and the potential ip with respect to the 
compact-open topology. 

Remark. Of course, Theorem|Sl has local variants 
but it is more convenient to formulate and prove it 
globally. 

Remark. Theorem|Sl looks correct for any smooth 
convex Hamiltonian h but the author has failed to find 
its proof in this case. 

7 Proof of TheoremUl 

Theorem|21 follows from Theorem|31 applied in the 
strip R'^ X [0, T] of the affine space-time to the poten- 
tials ip'^ and ip. Let 

X = (x,i), P=(p,a), Q = (q,T), 

h{p,a) = |p|V2 + a, 
e(x,<) = C7(r)(|xp+t2)/2 

where C{T) is the constant from Lemma^ According 
to Lemma^ the potentials ip'^ satisfy the conditions of 
Theorem|Slin the strip. Besides, 

where Vip'^ is the usual gradient of the smooth poten- 
tial ip'^. Theorem^] shows that the limit potential ip 
satisfies the conditions of Theorem|31 as well and the 
convergence ip'^ ^ (p as v ^ -1-0 is uniform in the strip. 
So, the derivative ip'^ ^ exists and it remains to show 
that 

V,,(^(x,i) = (u(x,t),l) (8) 

where u(x, t) is the center of the minimal ball contain- 
ing the set fcx.t- 

Remark. Formally speaking, we cannot apply The- 
orem|Slin the strip but we can always extend our poten- 
tials up to functions on the space-time satisfying the 
conditions of Theorem|31 

In order to show the equality (jSJ we can use the both 
forms of the definition /i-gradient from Sectional 

Hamiltonian form: According to jSJl, the set 2?x.t 
of the sub-differentials of the limit potential (p at the 
point (x, t) is the convex hull of the set 

{ (p, a) I p e K^^u \p\V2 + a + U{ic, t) = O} . 
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Hence, the Hamiltonian h{p,a) — |pp/2 + cr attains 
its minimum on Px,t at some point (u(x, i), cr*) where 
u(x, t) is the center of the minimal ball containing the 
set Kx^t- But substituting into © r = and compar- 
ing with the formula we get that the sets fcx,t and 
isTx^t has the same convex hull. Therefore, u(x, t) is 
the center of the minimal ball containing the set fcx,t 
as well, and we get 

V;,(^(X,0 =/iu(x,t),.. = (U(x,t),l). 

Lagrangian form: The Lagrangian 

r |q|V2 if r = 1 
^(^'^^ = 1 +00 if T^l 

is the Legendre transformation of the Hamiltonian 
h{p,a) = \p\'^/2 + a. This means that the Lagrangian 
form of the definition of the /i-gradient from Sectional 

is nothing but the principle 

8 Proof of TheoremEl 

According to Theory of ODE, if the potential ip is 
smooth then the Cauchy problem has a unique solution 
and the boundedness of the differentials of the potential 
ip guarantees that it is defined globally. 

Moreover, our Cauchy problem has a solution if the 
potential ip is homogeneous of the degree 1 and con- 
cave. (The last word means that the function —ip is 
convex.) This key observation is formulated in the fol- 
lowing lemma. 

Lemma 3. If p is a concave homogeneous function: 

piXX) = Xp{X) VA>0, 

then the trajectory X : [0,+oo) M™, X{t) = 
fS/hfiO) is a solution of the Cauchy problem 

X+{t)=VhpiXit)), X{0) = 0. 

This lemma is true for any smooth convex Hamil- 
tonian; in fact, the following proof uses its smoothness 
only, but the convexity is needed for the uniqueness of 
the ^.-gradient. 

Proof. Firstly, Vx C Vq for any point X € 
because our potential ip is concave and homogeneous 
of the degree 1. (Any sub-differential at any point is a 
sub-differential at 0.) 

Secondly, if the smooth Hamiltonian h attains its 
minimal value on Vq at a point Pq G Vq then Pq G 



T^xit) for any t > 0. It immediately follows from the 
equality 

ip{X{t)) = Po ■ X{t) 

because ip{X) < Pq ■ X . The last equality follows from 
(0. Independently on this reference: X{t) = th'p^ and 
for i > 

^(t'^'po) = min^ P ■ f^'po = Po ■ ih'p^r 

Indeed, for any P G we get {P—Po)-h'p^ > because 
Pq is a minimum point of the smooth Hamiltonian h 
on the convex set Vq. 

Therefore, the inclusions Pq G T^x{t) C Vq show that 
Pq is a minimum point of the Hamiltonian h on the set 
T^xit) which implies that \7h'f{'^{t)) — Vh¥'(0). □ 

Let $6 be the space of all potentials p with bounded 
sub-differentials such that the differences e — ip are con- 
vex functions (i.e., <I>e consists of the potentials satis- 
fying the conditions of the theorem) and 

g : a>e X M+ X M'" ^ M™, {p, t, Xq) ^ g^(Xo) 

be the mapping sending a potential ip, a time t > 0, 
and an initial point Xq to the value of the solution of 
our Cauchy problem at the time t. In other words, the 
trajectory 

Xit)=g'^iXQ) 

is a solution of the Cauchy problem. 

A priori, the mapping g can be many-valued and 
defined not everywhere. Let D{g) C $e x M+ x is 
the domain of definition of the mapping g. 

Lemma 4. // the Hamiltonian is the sum of a linear 
form and a positive semi-definite quadratic form: 

h{P) ^ h,{P) + h2{P) 

then the mapping g is one-valued and continuous on 
£'(g) with respect to the compact-open topology in the 
space $e- 

Remark. In the case h{P) = |Pp/2 the proof of 
Theoremn] proves, in fact, Lemma^Jas well. The key 
place is the inequality (|SJ|. 

Proof. We are showing the continuity of the map- 
ping g at a point ((ys,, t,, X,) G x R+ x M™. Let us 
consider an open bounded set K consisting the trajec- 
tory 

X,{t) ^ gl^iX,), tG[0,4] 

- it can always be done because the /i-gradient of the 
potential ip^ G $6 is bounded. Let an open bounded 
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convex set M contain the closure K of the set K and 
Lp be any potential such that 

\ip{X) - ip,{X)\ < e VXgM. (9) 

According to TheoreniQl the last condition guarantees 
that the sub-difFerentials of all such potentials at all 
points of K are bounded by a constant B. 
Let us consider another trajectory 

xit) = gl{x), xeK 

and introduce the following notation: 

SX{t) = X{t) - X,{t), SP{t) = Vit) - V,{t), 
where 

are the sub-differentials of the potentials ip and 93* 
where the Hamiltonian h attains its minimal values. 
The key inequality 

6r{t) ■ dX{t) < 2e + 2 e{6X{t)) 

is almost the inequality Like there we have 

ifi^iX) - ifi^X^) <V*-{X~ X^) + e{X - A",), 

if{X^) - ip{X) <V-iX,-X)+ e(X - X^). 
Adding these inequality and taking into account that 

-2e < ip^X) - Lp^{X^) + ip{X^) - ip{X) 

according to lO, we get our key inequality. 

Our further proof is coordinate. Let A" — (Ai,A2) £ 
= xR'"^ Xi e M"i, X2 e M"'^ be affine coor- 
dinates such that in the dual coordinates (Pi, P2) G T* 
the Hamiltonian has a canonical form: 

h{Pi, P2) = P ■ Q + \P2\V2, QeT. 

Then our trajectories satisfy the differential equations 

X+it)=Q + iO,P2{t)), X+{t) = Q + iO,P,2{t)) 

because h'p ^ Q ■ dP + P2 ■ dP2. Hence, 

SX+ (t) = 0, SX+ (t) = 5V2 (t) . 

Let \Xi - A,i| < Qfi, \X2 - A,2| < a2, 

R{t) ^\SX2{t)\^, 2e{SX) <C\SX\^ 

for some C > 0. Then 

\SXiit)\ = |Ai - A,i| < ai 



and the above key inequality gives that 
SV{t) ■ 5X{t) <2e + C 

Hence, 

R+{t)=2 5V2it)- 6X2(1) < 

< 4:e + 2C\6X{t)\^ - 25Vi{t) ■ 5Xi{t) < 

< Ae + 2C al + 2C R{t) + A B ai. 

Solving the differential inequality and taking into ac- 
count that R{0) < a^, we get: 

2Ct _ 1 

R{t) < {2e + Cal + 2Bai) + e^^K 

Therefore, as (e, ai,a2) +0 

\6Xi(t)\^0, R{t)~^0 

if t < . This proves the continuity when time is fixed: 
t ^ t^,. But if i < we can use the estimate 

\X4t)-X4U)\ < {\Q\ + B)\t-U\, 

and in the case t > t^, — 

\X{t)~X{U)\ < {\Q\+ B)\t~U\. 

□ 

Now, in order to prove Theorem|5l it is enough to 
show that £>(g) — x M+ x K™. Of course, we can 
uniquely define a continuous mapping 

: $e X M+ X M™ ^ K", 0b(g) = g. 

Indeed, D{g) is dense in <i>e x M+ x M™ because all 
smooth potentials form a dense subset in <i>e - that is 
shown with the help of the standard smoothing. 
It turns out that any trajectory 

X{t) = g*^(Ao) 

is a solution of our Cauchy problem and, therefore, 
g = g. In order to show that we have to check 

X+(0) = Ao and X+{t) = \7h{^) (X(t)) . 

Of course, the first equality follows from the obvious 
fact B° = id which is implied by g" = id. Why is the 
differential equation satisfied? 

The point is that the class of our Cauchy problems 
is invariant with respect to adding constants to poten- 
tials, translations of A, positive shifts of t, and simul- 
taneous dilations of the graphs of potentials, t, and A. 
Hence, the mappings g and q are invariant with respect 
to these transformations as well because the latter are 
continuous in the spaces M+, and M™. In the terms 
of the mappings g* : M™ M™ it means the following. 
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1. Adding constants to potentials: fl^+const = 0^- 

2. Translations of X: Q%(x-Xo) = Q%(x) + ^o- 

3. Positive shifts of t: = g^i o 

4. Simultaneous dilations of the graphs of potentials, 

and X: 

Q\%{x/\){^^o) = ^Q%(x){^a), A > 1. 

(The last inequality guarantees that X(p{X/\) 6 
$e if G $e.) 

The invariance property 3 implies that the differen- 
tial equation is enough to be checked only for t = 0. 
According to the invariance properties 1 and 2, we can 
suppose that Lp{0) — and Xq — 0. Besides, the in- 
variance property 4 implies that 

0W/A)(O)-AX(t/A), A>1. 

Hence, 

X+(0)= hm AX(1/A)=4,(0) 

because X(0) — Xq ^ and 

Xip{X/X) (p'{X) in $e as A ^ +oo 
in consequence of (p{0) — and TheoremsOlandE] But 

according to Lemma|31 that completes proving Theo- 
rem|S| 

9 Computations for Figures 

Figurc|21shows how particles move around nodes and 
end points of the shock. Figures0HI|demonstrate what 
can generically happen with a cluster moving in the 
shock. In order to get all of these Figures except the 
one containing an end point (on the right of Figure!^, 
we use the following procedure for fc = 3 or 4. 

Let in a neighborhood of some point (x»,t») the 
limit potential be presented as the minimum of k 
smooth solutions of the Hamilton-Jacobi equation: 

^(x, t) = min {ip^{x,t),. . . ,ip'' (x, t) } 

where (y9(x,, i,) — ip^{x^,t^) = ■ ■ ■ = i^'^(x,, t,) and 

ipl{x,t) + |V^''(x,t)|V2 + L/(x,i) = 

for all i — l,...,k. Therefore, the derivative (p'^ — 
Vx. t, of the limit potential at the point (x, , t,) is pre- 
sented by the formula 

fli^l^T) = min {pi • q - r|pip/2} - [/, r 



where Pi = V(/5'(x»,t,) and [/* = C/(x,,t*). Let us 
linearize the shock and the motion of particles around 
the point (x^,<^) considering the derivative ip'^ instead 
of the limit potential itself. The shock of the derivative 
and the velocities of the particles are described with the 
help of the last displayed formula and Theorem|21in the 
following way. 

Let any three of the momenta pi, . . . ,pk do not be- 
long to the same line and any four of them do not 
belong to the same circle. 

For T < 0. If the minimum of the values jp^r — q|, 
I = 1, . . . , fc is attained only for one index ii then the 
particle q is outside of the shock and its velocity is 
Pi^. If this minimum is attained for two indices ii,i2 
then the particle q is being at a smooth point of the 
shock and its velocity is the midpoint of the segment 
[pijjPij]. If the minimum is attained for three indices 

^2, ^3 then the particle q is being situated at a node 
of the shock and its velocity is the center of the minimal 
disk containing the momenta p^^ , p^^ , pi^ . 

For T > 0. If the maximum of the values \piT — q|, 
I = 1, . . . , fc is attained only for one index ii then the 
particle q is outside of the shock and its velocity is 
Pij. If this maximum attained for two indices ii,Z2 
then the particle q is being at a smooth point of the 
shock and its velocity is the midpoint of the segment 
[Pii I Pi2]- If the maximum is attained for three indices 
ii, *2, *3 then the particle q is being situated at a node 
of the shock and its velocity is the center of the minimal 
disk containing the momenta p^^ , pi^ , pi^ . Besides, the 
center of the minimal disk containing all of the points 
PiT, . . . , pfcT belongs to the shock and the unique tra- 
jectory coming from the origin (q, r) = 0. 

Remark. If r > 0, the shock of the derivative is the 
Voronoi diagram of the points pir, . . . , p^r. If t < 0, 
the shock is the set being analogous to the Voronoi dia- 
gram of the points Pit, . . . , p^r but defined by multiple 
maxima of the distance (not minima). 

Let fc = 3 and the momenta pi, p2, and ps do not 
belong to the same straight line. Then the shock of 
the derivative ^p'^, has a node at the center of the cir- 
cle containing the points rpi, Tp2, and rpa and the 
velocity of this node is the center of the circle passing 
through the momenta pi, P2, and P3. But the velocity 
of the particle situated at the node at a given time is 
the center of the minimal disk containing the momenta 
Pi, P2, and P3. Of course, these centers do not always 
coincide. 

Namely, there are two generic possibilities: the tri- 
angle with the vertices pi, p2, and pa can be ob- 
tuse or acute. If the triangle is obtuse then the ve- 
locities of the node and the particle are different and 
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the particle leaves the node. Otherwise, if the triangle 
is acute then the velocities coincide and the particle 
stays at the node. It is convenient to apply the above 
procedure in a frame of reference connected with the 
node. In this case the shock does not change with 



time, |pi 



IP2I 



IpsI , and we get the shock and 



velocities shown in Figure|31on the left (the triangle is 
obtuse) and in the middle (the triangle is acute). 

Remark. There is another simple explanation of 
the difference between these possibilities. Namely, in a 
frame of reference connected with the node the deriva- 
tive ip'^ has no extremum on the left and does have a 
maximum in the middle of the Figure|3| 

At separate times our triangle can become right - it 
is shown in Figures 0] and [S] on the left. 

Let A; = 4, and pi,...,p4 is a generic configura- 
tion of momenta. It means that they do not belong 
to the same circle and any three of them do not be- 
long to the same line and do not form a right triangle. 
Such generic configurations have many connected com- 
ponents, the components with four obtuse triangles are 
called totally obtuse. A component is called narrow if 
the boundary of the minimal disk containing all the 
momenta passes only through two of them. A com- 
ponent is called wide if the boundary of the minimal 
disk containing all the momenta passes through three 
of them. Of course, any totally obtuse configuration 
is narrow. Each component, except the totally obtuse 
ones, defines one of the types of behavior of clusters 
during the fifth or sixth transitions from Figure[21- all 
these types are shown in the middle and on the right 
of Figures EHZI Namely, taking a configuration from 
a connected component we apply the above procedure 
which shows the following. 

Each triangle formed by three of the momenta de- 
fines a node before or after the transition. If this trian- 
gle is acute then there is a growing cluster at the node. 
If the configuration is totally obtuse then there are no 
clusters at all and we ignore it. If the convex hull of 
the momenta is a triangle then the fifth transition oc- 
curs, if the convex hull is a quadrangle - the sixth one 
does. If the configuration is not totally obtuse then 
there is a cluster after the transition - its trajectory 
comes from the origin. If the configuration is narrow 
then this cluster is stable, if it is wide then the cluster 
is growing. 

It remains to compute the velocities around an end 
point shown on the right of Figure|3| Let 

^(x,t) =mm{F(C,x,i)} 



Hamilton-Jacobi equation: 

Ft{^, X, t) + I V,F(e, X, i)| 2/2 + C/(x, t) = 

such that the function F{-, x», i») has the simplest min- 
imum from the degenerate ones at the point ^ = 0. 
(In Singularity Theory this minimum is called A3.) It 
means: 

F{^,^,t) - F{0,^,t) ^ 

k k 

= Aa^ + 2'^B, a^h, + ^ C^Jb,b■i + 

1=1 ij = l 

k 

+ a(q, T)a + /3(q, T)a'^ -f ^ 7,(q, r) 6, + . . . , 

^ = (a,5i, . . . ,6fe), q = x-x,, t = t - t^; 

where the first part of the right side is a positive definite 
quadratic form of a^, 61, . . . , 6^; a, /3, 71, jk are 
linear forms of (q, t); and the dots denote the high 
order terms. This minimum representation gives us 
the following. 

1. In a neighborhood of the point (xa,,f*) the shock 
of the limit potential ip in the space-time is ap- 
proximated by the semi-hyperplane 



a(q,T) = 0, 



det 



//3(q,T) 7i(q,T) 

Bi Cii 



7fc(q,T)\ 

Cifc 



Ckk / 



< 0. 



V Bk Cfei 

2. If p, = VxF(0,x,,i,) then 

<<=x.,t.(q;'^) = p* • q - t|p*P/2 - u{-K^,t^)T. 

Hence, Theorem|2 implies that u(x»,i*) = p* and 
the space-time vector (p,, 1) is tangent to the tra- 
jectory starting at the point (x, ,t,). 

3. Differentiating the Hamilton-Jacobi equation we 
get: 



Fat + VxF • Vxi^a = 0, Fba + VxF ■ VxF. 



0, 



VxF • V^Fa. 



'VxFa • Vxi^a. 



After substituting ^ = 0, q = 0, and t = these 
equalities give: 



where -F is a family of smooth solutions of the 



a(p„l) = 0, 7<(P*,l)-0, /3(p„l)<0. 
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Therefore, the trajectory starting at the end point 
(x*,^*) of the shock goes inside it because its tangent 
vector (p*, 1) belongs to the semi-hyperplane from the 
first observation. In order to see that, we have to take 
into account the third observation and that 



according to the positive definiteness. 
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